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Correlation-function expressions are derived for the coherent nonUnear response of molecules to 
three resonant ultrafast pulses in the x-ray regime. The ability to create two-core-hole states with 
controlled attosecond timing in four-wave-mixing and pump probe techniques should open up new 
windows into the response of valence electrons, which are not available from incoherent x-ray Raman 
and fluorescence techniques. Closed expressions for the necessary four-point correlation functions 
are derived for the electron-boson model by using the second order cumulant expansion to describe 

Othe fluctuating potentials. The information obtained from multidimensional nonlinear techniques 
could be used to test and refine this model, and establish an anharmonic oscillator picture for 
, electronic excitations. 
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INTRODUCTION 



The development of bright attosecond soft and hard x-ray sources such as the free electron laser has triggered 
considerable interest in all-x-ray nonlinear spectroscopy |Q 0, 0, H- In resonant optical techniques in the visible, a 
photon is tuned to high frequency 2eV) electronic transitions, but a wealth of information is provided on lower- 
frequency (< 0.4eV) nuclear (vibrational and phonon) degrees of freedom which are accessible through multiphoton 
''O ' (6-g- Raman-type) resonances with differences (and higher combinations) of visible photons. In a completely analogous 
^ ' manner, combinations of x-ray photons resonant with high frequency ( keV) core transitions can probe the lower 
O frequency (< 50eV) valence electronic excitations. By exploiting this analogy, we can use the theoretical apparatus 
developed for probing electronic and vibrational coherences in nonlinear optics 0, to predict nonlinear x-ray signals 
and design new multiple pulse experiments. For example, the concepts underlying multidimensional techniques^ 
which provide extremely valuable information on optical excitations in molecular aggregates can be extended to probe 
correlations among multiple core hole states. 

Nozieres and De Dominicis had proposed a model Hamiltonian for resonant x-ray processes in metals 0. The ~ uj" 
threshold (Fermi edge singularity) behavior of x-ray absorption was calculated. The interplay of the large number 
, of electron hole pairs and the vanishing of the many-electron wavefunction overlap with and without the core hole 
^— I (Anderson's Catastrophe) can result in either a diverging {a < 0) or converging (a > 0) lineshape at threshold, as 
lO , shown by Mahan^, Jt, J^]. Signatures of this singularity in femtosecond optical pump probe spectroscopy of doped 
■ semiconductor nanostructures have been studied. Higher order radiative processes such as spontaneous emission 
(fluorescence or Raman) or nonlinear wave mixing provide more detailed insights through multipoint correlation 
functions ^5,, ,.12, 13, 14, 15] . The formal theory of x-ray response closely resembles its visible or infrared counterpart 
where the relevant correlation and response functions and possible time orderings have been studied extensively, 
It is also closely connected with the treatment of currents in open molecular systems coupled to electrodes. ^3 
In this paper we apply the density matrix Liouville space formalism^ to calculate time resolved x-ray four-wave 
Q ^ mixing signals and compare them with spontaneous emission and pump probe spectroscopy. Much current activit 
O ' is focused on the application of time resolved diffraction to probe structural changes such as surface melting, [l^ 
From a theoretical point these can be described using the existing formalism of diffraction by simply including the 
parametric time dependence of the electronic charge density; coherence does not play a role in these techniques. Other 
experiments have been carried out using a visible or an infrared pump followed by the absorption of an x-ray probe. 
These techniques provide x-ray snapshots of vibrational coherence or coherence among valence electronic states. [T9| 
Photoelectron spectroscopy provides additional novel ultrafast probes. [20II2H 

Pure resonant x-ray nonlinear optics of the type considered here can probe high frequency coherence of many- 
electron states involving core hole transitions and provide multidimensional real space pictures for the response of 
valence electrons to external perturbations. The necessary multipoint correlation functions may be calculated using 
several levels of theory, (i) Multiple summations over the many electron states of the valence system with N, N + 1, 
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and N + 2 electrons in the presence of zero, one, and two core holes respectively, (ii) The transition potential method, 
which uses a reference wstem with partially filled orbitals 1231, 153, l25l| (iii) many-body Green function perturbative 
techniques. |E llfl IIj, llS l2a l23 (iv) Replace the original Hamiltonian by the electron-boson model (EBM) and 
treat charge-density fluctuations as classical oscillations. |2a.l28L l29. 30, 31] 

Method (i) is conceptually the simplest but the most expensive numerically. The sum over state (SOS) expressions 
for nonlinear response functions allow us to employ any level of quantum chemistry for computing electronically 
excited states. TDDFT, for example, provides a relatively efficient way for computing a large number of electronically 
excited states. The relevant states are determined by the pulse bandwidth, e.g. 3.75eV for a lOas pulse. Method 
(ii) Represents systems with different numbers of core holes by varying the occupation numbers of a single set of 
reference orbitals. This approximate method works well for core level spectroscopies of small molecules, and may be 
extended to the nonlinear response. Method (iii) was used by Nozieres and coworkers to compute XANES and x-ray 
Raman spectra. 0,0] It is formally exact, and allows the development of powerful approximations. The EBM is by 
far the simplest to implement, since it is exactly solvable by the second order cumulant expansion. The model is not 
expected to apply for small molecules where the core electrons added to the valence band by the absorption of an 
x-ray photon need to be treated explicitly. It has been tested and found to work quite well for solids (semiconductors 
and metals). Plasmon satellites have been predicted in x-ray photoemission.p. llsL l3^ and inelastic scattering. [3^ 
They are less clearly seen in x-ray absorption, but have indirect signatures in e.g. the oscillator strengths. Using 
the EBM, electronic excitations are described as anharmonic oscillations whose parameters may be extracted from 
coherent nonlinear x-ray techniques. 



II. THE NONLINEAR X-RAY RESPONSE 

We start with the Mahan-Nozieres-De Dominicis(MND) HamiltonianfTl H I2 II33L Issf 

-ff = ^ Cfecj^Cfc -I- ^ e„dJid„ + ^ Vkk',nclck'dnd'l (1) 



kk' 



where , c?Jj are the Fermi creation operators for a valence and a core electron respectively, e„ is the energy of core 
state 71, and Vkk'.n is the potential acting on the valence electrons due to the nth core hole. 
The dipole interaction with the x-ray field in the rotating wave approximation (RWA) is 

= t^^'"' *)-^« + ^* ^)^"] (2) 

n 

where 

BI^'Y M/cnC^dn, -S" X! f^nkCkdi (3) 

k k 

are creation (annihilation) operators for core-hole excitons. fikn is the dipole matrix clement between the n'th core 
orbital and the fc'th valence orbital, and E{r, t) is the complex field envelope (see Eq. (|25|) '). The possible transitions 
between the zero, one and two core hole states are shown in Fig. 1. 

The j'th order molecular response to the x-ray field is described by the induced polarization 0: 

P(^)(r,i) = rr(Pp(^)(t)) (4) 

where Tr(. . .) denotes the trace and 

P = ^(S„ + Bi) (5) 

n 

is the polarization operator, p'^^^ (t) is the density matrix describing the state of the molecule obtained by solving the 
Liouville equation to j'th order in the field. 

f^-^[H,pl (6) 
where the total Hamiltonian H = Hq + Hint , is the material Hamiltonian. 
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Eq. (0J can be expanded as 

p(-')(r,<) = / dt,J dt,,_i... I dti5(-')(t„,t„_i,...,ii) X 



^0 

E{r, t ~ tn)E{r, t - - tn-i) . . . E{r, t - t„ - t„„i ... - h). (7) 

where S^-'-' in_i, . . . ,ti) denotes the j-th order response function and t„ = t„+i — t„ is the time delay between 
two consecutive interactions with the x-ray field (see Fig. 2). 

The first order polarization is related to the linear response function S"*-^^: 

POO 

P(i'(r,<)= / dhS^'\h)E{r,t-h), (8) 



S^^^ is a second rank tensor with respect to the polarization direction. For clarity we do not use tensor notation. 
Here 

s^'Hh) = t9{h)[j{h)-r{h)], (9) 

^(''i) = (^('•i)^(O)) is 9- two point correlation function and the polarization operator is given in the Heisenberg 
representation: 

P{h) - eM^Hoh)PeM-iHoti), (10) 

6{t) is the Heavyside function {9{t) = for t < 0, d{t) = 1 for t > 0) which represents causality, and the angular 
brackets (• • • ) denote the trace over the ground state density matrix: 
The third order polarization is given by 

P(3)(r,t) = JJJ dt3dt2dhS^^\h,t2,h) X 

E{r, t - h)E{v, t~h- h)E{r, t-h-t2~h). (11) 

The third order response function is similarly given by a sum of eight terms, each representing a distinct Liouville 
space pathway j^: 



s{h, h) = i^e{t3)e{t2)e{h) J2 [Rpih, t2, h) - R^hMM)] , (12) 

where 



RiihMM) = F{h,h + hM+t2 + h,{)), 

R2it3,t2,h) = F{0,h+t2,h+t2+t3,h), 
R3it3,t2,h) = F{0,h,ti+t2 + t3,h+t2), 

R4{t3,t2,tl) = F{ti+t2+t3,ti+t2,h,0), (13) 

and the four-point correlation function is given by: 

F{n, T3, T2, ri) = {P{Ti)P{T3)P{T2)Pin)). (14) 

The polarization may be alternatively expressed in the frequency domain using the susceptibility tensors: 

X^^^(-Wa;Wa) = X^^H^a) and X^^H-^s;UJa,UJb,UJc)- 

nOO 

pW(r,i)=/ dLJaexp{-iLJat)x^^H^a)£{r,iOa). (15) 



P^^\r,t) — J JJJ du!sdLUadu!bdu!cexp{~iujst) X 

xx^^\-uJs;uja, uJb, ujc)S{r, uja)S{r, ujb)£{r, cUc)- (16) 
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S{r,uja) is the x-ray field in the frequency domain: 

/oc 
dTE{r,T)exp{iujT). (17) 

The response functions and the nonhnear susceptibihties are related by a Fourier transform: 

X^^\uja)= dhS^^\h)exp{iiJati), (18) 



exp{i{uJa +UJb+ ljJc)h + i{ljJa + i^b)t2 + iuJatl), 

where uj^ = uja + i^b + and the sum runs over all 3! = 6 permutations of cJq, o^h, Wc- 

III. EXPANSION IN MANY-ELECTRON EIGENSTATES 

Expressing the two-point correlation function Eq. (PJ using Eq. (|SJ) we get 

J(i)=5](B„(i)Bt(0)) (20) 

n 

Similarly, upon the substitution of Eq. ((SJ in Eq. (|14f) we find three contributions F = Fi + F2 + -F3 where 

i^l(T4,T3,T2,ri) = ^(B„,(T4)Bt,(T3)B„(T2)Bt(Tl)), 

n.'m 

i^2(r4,T3,T2,ri) = ^ (B„(T4)B„(T3)Bt,(T2)Bt (ri)}, 

F:i{T^,T3,T2,T,) = J2iB,n{n)BM3)Bl{T2)Bi{n)). (21) 

Fi only contains transitions to and from the ground state g—^n^g—fm^g^ and depends on either one (n = m) 
or two (n 7^ m) core holes. Its sensitivity to correlations between core holes stems from an interference between two 
pathways that lead to the same valence electron-hole pair via the two possible intermediate channels (core hole on n or 
m). At no point along the path do we have a state with two core- holes existing simultaneously. F2 and F3, in contrast, 
are intrinsically cooperative since they also include transitions among the excited states g n nm —>■ m g and 
depend on two-core-hole (two exciton) states. This may best be seen using the Liouville space pathways (Eq. (|13|) ) 
displayed in Fig. |3| 

The evaluation of these matrix elements requires the many-electron wavefunctions {ip^) of the original molecule 
with N valence electrons, where v = g is the ground state and v = e, / . . . are valence excited states. In addition 
we need the valence N + 1 electron wavefunctions calculated in the presence of the nth core hole, \ipnt'^) ^.nd iV -t- 2 
electron wavefunctions calculated in the presence of two core holes at n and m, lipnmt)- The corresponding energies 
will be denoted Ej^ , + ^ -I- il„ and E^^l + ^ln + ^m respectively. Here f2„ is the core-hole excitation energy whereas 
E^ is the energy associated with valence electrons. Expanding in these states, the four point correlation functions 
(Eqn. 1)21(1 ') assume the form 

Fi(r4,r3,r2,ri) ^ J2 « 1 5™ (^4) |V'^+3')(V'™j:3 1^^.(^3)10 

X «|S„(r2)Kii)«+/|i?t(ri)|Vf) 

F2(T4,T3,r2,n) = E E (V-f |S«(^4)Ki')«i'|5™(^3)|V'^+l,) 

X i^nr^^X |S,^«(^2)|V^^i')«i'|Si(ri)K) 
F3{r,,Ts,T2,n) (V^f |S™(^4)|V'^+3')(V'™+.3l^n(^3)|V'lr.l.) 

X (V-^r.l. |fiL(r2)|V^i^)(^^+/|Si(ri)|0 

(22) 
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These can be expressed as 

^^1(7-4, T3, T2, Tl) = ^ exp|^ - iQrnT43 - i^nT21^ ^ CXp - iE^'^lr^s - iE^^Tz2 - iE'^ ^,^^T2l 
^2(1-4, T3, T2, Tl) = ^ exp(^ - in,nT32 - i^,iTii^ ^ exp - iE^+^^Ti3 - iE^^^^^T32 - iE^^^T2l 



n^m 1^1 1^2 



^E,n,^:,ri3 - lE,^^^^^T32 - lE^ ^^ T2I 



7N+2 



where Tij = — and i,j = 1, 2, 3, 4. For the hnear response we have 



J{r) = ^ |«rii?ilOrexp(-zl]„r - z<+V) 



(23) 



(24) 



A finite core-hole lifetime can be added by setting 



2 '~im 1 ^71 



nr. 



j7„. 7 provides a time window 



for the experiment. Typically it is ~ 0.375eV which corresponds to ^ 10 fsec window. Only higher frequencies and 
faster processes than this window can be probed by resonant x-ray techniques. Information about multiple core-hole 
dynamics can be also extracted from frequency domain x-ray four wave mixing. [3^ 



IV. COHERENT MULTIDIMENSIONAL SIGNALS 



We consider a sequence of x-ray pulses (Fig. whose electric field is given by: 

4 

E{r, i) = ^ Ej (t) exp {ikjr - iuj^t) + c.c. (25) 
i=i 

Here Ej (t) is the slowly-varying complex envelope function of pulse j with carrier frequency coj and wavevector . 
c.c. denotes the complex conjugate. Most generally, a third-order process requires four external fields: three (j=l, 2, 
3) interact with the system and the fourth, heterodyne, field {j—4:) is used for the detection. 
To calculate the signals we expand the nonlinear polarization in k space: 

P(3)(r,i) =J2p(^^t)eM^sr), (26) 

where the possible wavevectors are k^ = ±ki ± k2 ± ka . 

We shall consider well-separated pulses where pulse 1 comes first, followed by 2 and finally 3|2lThree signals are 
possible for our model, k/ = — ki + k2 -I- ka , k// = ki — k2 -I- ka , and k/77 = ki -|- k2 — ks . The polarizations responsible 
for the these signals, obtained by invoking the RWA (i.e. neglecting highly oscillatory off- resonant terms), are given 

by 
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Fi{ti,T2,T4,T3) + Fi{ti,T3, T4,T2) - F2 {ti , T4 , T3 , T2) - (ti , T4 , T3 , T2) 
X EliTi)E2iT2)E3iT3) 




Fi{t2,T3,T4,Ti) + Fi{t4,T3,T2,Ti) - F2 (t2 , T4 , Ta , Ti ) - F3 {t2 , T4 , T3 , Ti) 
X ^l(ri)^2*(^2)^3(T3) 




F2{t4,T3,T2,Ti) + F3{t4,T3,T2,Ti) - F3 (ts , T4 , r2 , Tl ) - i^2(T-3,T4,r2,ri) 
X ^l(ri)^2(T2)S3*(^3) 

(27) 

For very short (impulsive) pulses, we can eliminate the time integrations and simply set ri ~ t — t3 — t2 — ti, 

(3) • 

T2 = t ~ t3 — t2, T3 ^ t — t3, and T4 = t. Pj will then depend parametrically on the three time delays ii, t2, and ia 

(Fig.m. 

The physical processes underlying each of these signals can be understood by using the Feynman diagrams shown in 
Fig. (PI 5] which depict the evolution of the valence electronic density matrix in the course of the nonlinear process. 
The two vertical lines in the diagram represent the ket and the bra (time goes from the bottom to the top), while 
arrows represent interactions with the laser pulses. A coherencedg) (m| or \m){g\) is created by the first pulse. The 
second pulse takes the system either to a single exciton population (|n)(n|), coherence (|n)(TO|), or to a two exciton 
coherence {\nm)(g\). The population and the coherence evolutions can then be probed by holding the second delay 
time, t2 fixed. The third pulse creates coherences either between the ground and one-exciton states or between one- 
and two-exciton states. The four diagrams contributing to k/ are shown in the left column. In all diagrams the 
density matrix represents a single-quantum coherence \g){n\ between the ground state and the singly excited state 
during ti. During ^3 it is either in the conjugate coherence \m){g\ (a and b), or in a coherence between the one and 
two exciton manifolds \mn){n\ (c and d). k/// is similarly described by the four diagrams ((i),(j),(k) and (1)) and 
shows double-quantum coherences between ground state and the two-exciton band \nm){n\ during the ^2 interval. 
During ^3 it has a single quantum coherence \n){g\ (i and j) and \nm){n\ (k and 1). k/, known as the photon echo 
technique, can improve the resolution by eliminating certain types of inhomogcncous broadening, k/// carries direct 
information regarding the coherence between the two exciton states and the ground state (double quantum coherence) 
so that its spectral bandwidth is doubled. [5l.l3^ 

Note that the absolute time arguments of F in Eq. (|27|l are not time ordered. However they are ordered on the 
Keldysh loop shown in Fig. Q in the following sense: for each diagram we can start in the bottom right, move up 
on the bra line and then down on the ket line. This will give the order of the r arguments for each of the 12 terms in 
Eq. (HZI). 

Within the slowly varying amplitude approximation the signal field is proportional to the polarization, Es{t) oc 
ipj-"'' (t) . The simplest detection measures the time integrated signal field intensity, and the third order signal in 
the kg direction is given by: 

/+00 
\P^^\t)\^dt. (28) 
-00 

This is known as the homodyne detection mode. Additional time resolution may be achieved by time gating which 

yields the absolute value of polarization itself Ihom{ii, ^2,^3) — \Ps'^\i)\'^ ■ 

In heterodyne detection the generated field Es{t) is mixed with a fourth field, E4{t) which has the same wavevector 
and the heretodyne signal is given by: 

/+00 
El{t)Pi''\t)dt. (29) 
-00 
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The time resolution is now determined by the heterodyne field, and the signal depends linearly rather than quadrati- 
cally on P, {t). By choosing different phases of the heterodyne field it is possible to measure separately the real and 
the imaginary parts of the polarization. 

A mixed time/frequency representation of the signal may be useful to reveal correlations in the system. For example, 
one can display two dimensional uji jio^ correlation plots for a fixed ti . 

//iet(t^i,i2,W3) = / / dtidt3lhet{ti,t2,t3)eyip{iujiti + iuJsts). (30) 

JQ JQ 

Time and frequency resolved signals and fields may be displayed using the Wigner spectrogram j^i 

/ + 00 
E; {t - t/2)Es (t + t/2) exp(ia;r)fir. (31) 

The spectrogram directly shows what fraction of the field energy is contained in a given time and frequency window. 
Integrating over the frequencies gives the instantaneous field energy 

Ws{t,uj)dijj ^2'K\Es{t)\^ (32) 

-oo 

while integrating over the time gives the energy density spectrum 



Ws{t,uj)dt=\£s{^)\. (33) 

The one dimensional projections of the spectrogram (Eqs. (|32|l and H33I) ) are known as marginals. 

To express the heterodyne signal, Eqn. (|29l) . in the Wigner representation we assume that the heterodyne field is a 
replica of one of the incoming fields in a nonlinear experiment, and expand the polarization to first order in this field: 

Ps{t) - / dTS'~^Ht,T)E^{T). (34) 



5 is a non-equilibrium correlation function of the system driven by all other fields. Defining the mixed time-frequency 
response function 

r + oo 

S^^^t,uj)^ S^^^t + T/2,t-T/2)exp{iujT)dT, (35) 



the heterodyne signal assumes the form [33, |3 



(• + 00 /• + 00 7 

IhetihM.h)^ dt —W4{t,oj)S(^^\t,uj), (36) 



2tt 



where W4(t, w) is the heterodyne spectrogram, (Eq. H31|ll. 

Eq. (|3(j|l is exact and holds for arbitrary field envelopes. For impulsive (very short) pulses the Wigner distribution 
is narrowly peaked at the time of heterodyne field and Eq. (|36|l reduces to 

hetitiMM) = S^^^\fi,fi) cc Im{El{fi)PM)}- (37) 

In the other extreme of ideal frequency domain experiments the spectrogram is narrowly peaked around its carrier 
frequency Cu^ and 

IhetihMM) = ^^'^^4,^4) (X Im[£l{uji)P,{uj^)\, (38) 

where 

/ + 00 
dTidT2S''^\Ti,T2) exp(ia;iTi + iuj2T2)- (39) 
-C30 
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V. PUMP-PROBE SPECTROSCOPY 

Pump-probe is the simplest third order technique and only requires two pulses. The signal defined as the difference 
absorption of the probe with and without the pump is related to the polarization at = ki + k2 — ki originating from 
two interactions with the pump (wi,ki) and one with the probe (w2,k2). The probe serves as the heterodyne field 
since the signal is measured in the probe direction. This technique may thus be viewed as self-heterodyne detection. 
This is an incoherent technique whereby the contributions of different molecules to the signal itself (rather than to 
its amplitude) are additive. 

We consider a sequential pump probe signal induced by short, well-separated pulses where the pump comes first, 
followed by the probe with a delay time of r. The signal is obtained from Eqs. (|27|l and 1)29(1 by combining the k/ 
and kji polarizations. We get Spp — S^^ -\- Spp, where 



/oc poo poc poo 

dt dh dt2 dtiE*(t-T + h)E2{t-T)El{t-t2)E^{t-t2-h) 

-oo Jo Jo Jo 

X exp[iuj2t3 + iLUiti] 

X Fi{t-t2-ti,t-T + t3,t-T,t~t2) 
+ Fi{t-t2-ti,t-t2,t-T,t~T + t3) 

/OC pOO pOC pOC 

dt dt3 dt2 dtiE;{t~T + t3)E2{t-T)Ei{t-t2)El{t-t2-ti) 

-oo Jo Jo Jo 

X exp[iLL!2t3 — iLUiti] 

X Fi{t~t2,t-T + t3,t-T,t-t2-ti) 
+ Fl{t~T,t~T + t3,t-t2,t-t2-tl 



(40) 



SE,{uji,U2;t) 



/oo poo POC POO 

dt / dt3 / dt2 / dtiexp[—iuj2t3 + iwiti] 
-oo Jo Jo Jo 

X E2{t + t3- T)El{t - T)E*{t - t2)Ei{t -t2- ti) 
X F2{t-t2-ti,t-T + t3,t-T,t-t2) 

+ F3{t-t2-ti,t-T + t3,t-T,t~t2) 

/OO pQO pOC pOC 

dt dt3 dt2 dtiE2{t + t3-T)E;{t-T)Ei{t-t2)El{t-t2-ti] 

oo Jo Jo Jo 

X exp[— 16^2^3 — iiuiti] 

X F2{t-t2,t-T + t3,t-T,t-t2-ti) 
+ F3{t ~ t2,t ~ T + t3,t - T,t - t2 - ti 



(41) 



These terms can be separated into two negative contributions, ground state bleaching (GB) and stimulated emission 
(SE) and a positive path, excited state absorption (ESA). The corresponding Feynman diagrams are shown in Fig. |S1 

Using the Wigner representation (Eq. 136(1 ). the pump probe signal can be expressed as an overlap integral of 
three functions: the pum p s pectrogram Wi{t',u!'), the third order response function S^^\t" ,uj" ,t' ,uj') and the probe 
spectrogram VF2(t",w") [3ll39l|: 



Ipp{ui,Ti,U2,T2,) 



dt'dt"duj'duj"W2{t\ w")S^^^ {t",u;", t',Lu')Wi {t',io'). 



(42) 



VI. FLUORESCENCE AND RAMAN SPECTROSCOPY 



Resonant x-ray emission is widely used in the study of core- hole transitions Q> We consider 

a molecule driven by an x-ray field with a complex envelop E{t) and carrier frequency lol- We shall write the 
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molecule-field interaction within the RWA in the form 

Hint = E(t)B^cxp (-iLULt) + E*(t)Bcxp (iujLt) . 



(43) 



B'^ {B) are the exciton creation (annihilation) operators and the dipole operator is given by /i —B + i?^. The time 
and frequency resolved flourescence spectrum is given by |4l) 



/oo /"t rt-\-T 

dT dri dT2{B{T2)B\t + T)B{t)B^ [ti]) 



(44) 



We can divide the response into time-ordered contributions which separate the Liouville space pathways into Raman 
and fluorescence typespj- The three pathways with r > and r < are complex conjugates. This gives 

/•oo pOO /"OO 

SF{cJL,i^s,t) = 2Re / dt3 / dt2 / dti 
Jo Jo Jo 

E{t ~ ti ^ t2 - t3)E*{t - t2 - t3)exp{iuJLtl + iuJst3)Fi{ti,ti + t2,ti + t2 + t3,0) 

+ E*{t^ti-t2~ t3)E{t -t2- t3)exp{~iujLti + iujst3)Fi{0, ti + ta, ti + ^2 + ^3, ti) 
+ E*{t - ii - ^2 - h)E{t - t3)exp(-iwi,ti - i{ujL ~ ^5)^2 + i^sh) 

X i^l(0,ti,ti+i2+<3,tl+<2)] (45) 

The three terms in Eqn. (|45|l come from i?i, R2 and R3 respectively of Eqn. IjlSI) . It is interesting to note that Sp 
only depends on Fi ; four wave mixing signals also depend on F2 and F3 , and therefore explore new regimes of Fock 
space not accessible by fluorescence. 

We now insert the complete basis of many-body states in Eq. H22|) . We denote the transition frequency between 
states \a) and \h) as LOab = — and Tab — \{la + lb) + Tafc is the corresponding dephasing rate. Here 7a is the 
inverse lifetime of state \a) and Tab is the pure dephasing rate resulting from frequency fluctuations. Assuming a c.w. 
field {E{t) = 1), we can carry out the time integrations and obtain 



SF{0Js,i^L,t) = ^'^^'-^ ( SliuJs,UJL) + Sll{uis,i^L) + SlIl{uJs , UJl) ) 



(46) 



where 



'5/(^5, Wi) 



(^,^+^|i3tJ^,^) «|i3„K+;) K) 
(wi/1 1/2 - + i^uiu^) {^1^1 1'3 + ^r^i ^3 ) (w^i g - -I- zE^j g ) 



(47) 



SlliuJS,^L) 



,.y y (V'^+^|j3t,|0 «|i?nK+;) K) 

'■^^^"^^3 (wi'ii'a -(^s + ir^i^j(w^i^3 + iF^,^3)(tJL - tj^gg + ir^3g) 



(48) 



Sin{uJs,ujL) 



E 



(C+^|i?J;JO «+/|i3t K) 



In the absence of dephasing, Tg^, = 0, (i^ = z^i, 1^2, t'a), these terms can be combined to yield 



27r 



E 



EE 



(49) 



(50) 



where 77 is an infinitesimal positive number. This is the standard expression for spontaneous emission spectral^. 
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VII. THE ELECTRON-BOSON MODEL 



So far we derived exact expressions for the multipoint correlation functions in terms of the many-electron states. 
A much simpler description can be obtained by replacing the valence excitations by a boson bath described by the 
operators , a| and adding the electron-boson model Hamiltonian [s^, |^ IssL 



Ck- 



(51) 



The first two terms represent the reference Hamiltonian for the noninteracting valence and core electrons. The last 
term is the potential induced by the fc'th core hole which causes a linear displacement of the bath modes. 



Uk 



n, m 



(52) 



Each time we act with B'' we add a valence electron. If the valence system is very large (e.g. the electron gas or 
a large metal nanoparticle) we can ignore the effect of the additional electron and only consider the added core hole 
potential. In this case we can set = c^, -B = c and the polarization operator becomes 



We then have 



P = ^(c, + 4). 

k 



B\{.t) = exp(zi/„T)exp(-ii7oT) 



(53) 



exp_ 



dr[/„(T) 



(54) 



where Hq consists of the first two terms in Eqn. (|51|l so that J7„ = i?„ — Hq is the fluctuating potential caused by the 
nth core hole. 



Similarly 



-B„(t) = exp^ 



-I I drUnir) 





(55) 



(56) 



exp_|_(-_j refer to positive (negative) time-ordered exponentials. The multipoint correlation functions may then be 
calculated using the second order cumulant expansion. ^3 

This is formally analogous to a multilevel molecular aggregate Hamiltonian with diagonal energy fluctuations whose 
nonlinear response was calculated in |43| . Assuming that the transition dipole operator does not depend on the 

bath coordinates, we obtain 



J(t) = (cxp(-i j dnUnin)'^ 



(57) 



J{t) = ^ |/Xg„pexp[-irj„T - ^5„„(t)], 



(58) 



flag is the transition dipole moment between the ground state and state a and fiat is the transition dipole moment 
between excited states a and b. The line broadening function g„„i is associated with energy level fluctuations: 



9ab{T) = \ ( dn f ' dT2 [{Ua{Tl)Ub{T2)) + {Ub{Ti)U a{T2))] 

^ Jo Jo 

The four-point correlation functions may be calculated by starting with j4,3| . 



(59) 
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Fi{t4,T3,T2,Ti) = ^^CXp+^-i J [/,„ (T)dT^ CXp_ J U„i{T)dT^ 

X exp_|_ ^— i J C/„(r)(ir^ cxp_ J Un{T)dT 
-p2(-r4,T3,T2,ri) = ^ ^exp_^ ^-i j [/„(r)dT^ exp+ ^-i ^ 



Ujn{T)dT 



X exp_ J U,n{T)dT^ exp_ y Un(T)dT 



X exp_ y C/,„(T)dT^ cxp_ y C/„(T)dT^ ^ 
Performing the cumulant expansion to second order . [4^ l47j we get 



exp[-«0,„T43 - iil„r2i - /i(t4, ts, T2,ti)] 



with 



with 



and 



with 



(60) 



(61) 



/l (n , T2 , T3 , r4) = gnn (t21 ) + (T43) + .9nm(T32) + Qnrn (r4i)- 

-ffnm (T3I )- 5nm(T42)- (62) 



-F2(T4,r3,T2,ri) = ^ 

f-^gnf-^gmf-^mg f^ng ^ 

exp[-irimr32 - iri„r4i - /2(t4, ts, T2, n)], (63) 



/2(n, T"2, T3, T4) = gmi('''4l) + 9mm{T32) — 5nm(T'2l) 

nm 9 ran 9 ran (T42). (64) 



^3(7"4,T3,T2,ri) ^ 

f-^gmf-^gnf-^mg f-^ag ^ 

exp[-zrimr42 - iri„T3i - /3(t4, rs, T2, ti)], (65) 



/3(t4, T3, r2,Ti) = g„„(T3i) + 5mm (T42) + 5„m(r32) - 5„m(r2l) + g„m(T4l) " gnm(T43)- (66) 

The fluctuating potentials enter the response through the spectral densities: 

1 /■+°° 

C"ai{uj) = -- diexp(zL^i)([f/,(0,C/fa(0)]), (67) 

where the expectation value and the time evolution are taken with respect to the bath Hamiltonian. These contain 
all relevant information about the fluctuations necessary for computing the nonlinear response of the system. gab{i) 
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are directly related to the diagonal spectral densities of the bath C" ab{^) 0: 



duj 1 — cos{ujt) 



coth 



+i 



2tt w2 
°° duj sin(ct;i) — tot 



2n 



\2kBT 



These spectral densities can be obtained from photooelectron spectroscopy or inelastic x-ray scattering, 
model for the bath is given by the overdamped Brownian oscillator spectral density: 



C"ab(l^) — 2Xab- 



^Kb 



^^ab 



(68) 
A simple 

(69) 



Aaa represents the magnitude of fluctuations of the energy of state a, while Xab represents the fluctuation of the 
coupling between states a and h. It can be observed in fluorescence as the time dependent Stokes shift. We further 
define the linewidth parameter A^j, = 2kBTXab- Substituting the overdamped Brownian oscillator spectral density 
we get in the high temperature limit: 



labit) 



^2( 



2TKb _ ^ Aofc (g^p(_A^^|^|) ^ ^^^i _ 



V A2 

\ a 



(70) 



Two dimensionless parameters, 77 and k, can be used to characterize the model and classify different regimes 
of energy fluctuations. The first, 77, defined by = Vab^aa^bb, represents the correlation of fluctuation 

amplitudes. —1 < r)ab < 1- These may be anti-correlated {rjab — —1), uncorrelated {rjab = 0) and fully correlated 
iVab = 1)- The second parameter, Kab = ^ab/^ab, is the ratio of the inverse time-scale of the bath to the amplitude of 
the fluctuations. It controls the lineshape; in the slow bath limit {Kab < 1) it has a Gaussian profile which gradually 
turns into a Lorentzian as Kab is increased |^ |48| . 



VIII. DISCUSSION 



Nonlinear core-hole spectroscopies could provide critical tests for the limitations of the electron-boson model. Even 
when the core holes are localized, the valence orbitals are usually delocalized and the third order response (Eqs. 
(|22J)) requires four sets of valence orbitals corresponding to l^i^"*""^), \ipn~^^) and lipnm^)- New insights could 

be provided on electron dynamics in conjugated molecules. |l4l Isq and mixed valence compounds. (49i| In the EBM, 
the effect of both core holes on the boson system is additive, (Eqn. t^52\i ). This is why the response only depends 
on two collective coordinates (fluctuation potentials), Un and Um- It is possible to extend this model in various 
ways. For example, if we add an extra coupling to Eqn. I|51|) H' ~ X^nm Kim,s(as + o,l)Cncl^CnCm, the dynamics 
will depend on a third collective coordinate, [/„„i = X^s ^im,s(as + ^l)- The bath displacement when both n and 
m holes exist simultaneously will then be non-additive, [/„ + [/,„ + [/„,„. The nonlinear response for this model 
may be calculated using the general expressions for a multilevel system given in ref ( |42j| ) . Fluorescence will not 
depend on Unm-, since it does not affect Fi. The valence electrons generally behave as anharmonic oscillators and 
the Hamiltonian may contain higher order terms in and a| . Coherent nonlinear x-ray techniques may provide new 
information about these anharmonicities in the same way that infrared multidimensional techniques probe vibrational 
potential surfaces. 6] The signal generally depends on many pulse parameters and can be displayed by various types 
of multidimensional correlation plots. For example, displaying the signals as a function of the time delays ii, t2 and t^, 
provides a direct look at valence electron wavepackets. In the electron gas these probe charge density fluctuations. The 
third order techniques considered in this article offer numerous new probes for electron dynamics. Complementary 
information is provided by second order techniques such as sum and difference frequency generation. 51] These can 
also be analyzed using the formalism presented in the article. IH^ 

Nonlinear core hole spectroscopy provides new ways for probing the dynamics and response of electronic valence 
excitations by controlled attosecond switching of external potentials. Femtosecond pulses introduced in the eighties 
allowed real time probing of nuclear motions. Attosecond x-ray pulses make it possible to watch electronic motions in 
real time. These techniques could help develop more realistic anharmonic Hamiltonians for the core- valence couplings, 
and test the validity of approximate model Hamiltonians. 
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FIG. 1: The core exciton states which enter the third order response. \g) is the ground state, |n) and |m) are single core-hole 
states, and \nm} is a double core hole state. 
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FIG. 2: (Upper panel) Laser pulse sequence in a three-pulse, four- wave mixing experiment. Three pulses ki, k2 and ks create 
the nonlinear polarization in the sample, which generates the now x-ray field in the direction kg = ±ki ± k2 ± ks. k4 = k^ is 
the heterodyne field. (Lower panel) The peak ordering and time intervals. 
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FIG. 3: Double-sided Feynman diagams representing the Liouville-space pathways contributing to the third-order response in 
the rotating wave approximation. The level scheme is given in Fig. Shown are the diagrams contributing to the four-wave 
mixing signal generated along the various possible directions: k/ — — ki -I- k2 -I- ks, k// — ki — k2 -f ks, and k/j/ = ki -I- k2 — ks. 
Diagrams (a), (b), (e), and (f) correspond to Fi and only include one-exciton states. All other diagrams also involve two-exciton 
states, (i), (c), (g), and (k) correspond to F2 and (j), (d), (h) and (1) represent ^3. 



18 



FIG. 4: Each of the 12 non-time-ordered correlation function expressions for Fj (Eq. (|2ZJ) may be obtained from the corre- 
sponding Feynman diagram (Fig. O using time ordering on a Keldysh loop. This is illustrated here for k/(a). (Fi(ri, r2, r4, rs)) 
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FIG. 5: Feynman diagrams representing the eight contributions to the sequential pump-probe spectrum. Diagrams (l)-(8) 
correspond respectively to the eight terms in Eqs. MUll and 14111 . (l)-(4) represent the F\ contributions (Eq. 14Ull ). which can 
be divided into ground state bleaching (GB) and stimulated emission (SE) type. (5)-(8) are the F2 and contributions (Eq. 
I4II 1'). which represent excited state absorption (ESA). 



